S-wave and p-wave scattering in a cold gas of Na and Rb atoms 
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Using improved experimentally based X 1 E + and a 3 E H 



the variable phase method to compute new data for low energy scattering of 23 Na atoms by 8E 



molecular potentials of NaRb, we apply 

"Rb 

atoms and s 'Rb atoms. These are the scattering lengths and volumes, numbers of bound states and 
effective ranges, which we use to obtain the low energy spin-change cross section as functions of the 
system temperature and the isotope masses. From an analysis of the contributions of s-wave and 
p-wave scatterings to the elastic cross section we estimate temperatures below which only s-wave 
scattering is dominant. We compare our quantal results to data obtained from the semiclassical 
approximation. We supply evidence for the existence of a near zero energy p-wave bound state 
supported by the singlet molecular potential. 

PACS numbers: 03.65.Nk, 34.10.+X, 34.20.Cf 
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I. INTRODUCTION 

The properties of ultracold trapped ensembles of atoms 
are governed by low-energy collisions that occur within 
such systems Of particular interest is the study of 
heteronuclear alkali-metal dimers because it is now pos- 
sible to generate dual species Bose-Einstein condensates 
(BEC) and produce Bose-Fermi ensembles from such bi- 
nary mixtures 0, 0, 0, HI • At sufficiently low tempera- 
tures a few scattering parameters suffice to characterise 
collisions between two atoms of different species. These 
are the s-wave scattering length, a s , the p-wave scattering 
volume, a p , and the effective range, i? ff, for scattering 
via the singlet and triplet ground states molecular poten- 
tials; but their accurate computation remains non-trivial. 

It is important to obtain precise scattering data since 
they give valuable information on the processes that al- 
low or not the production of a stable condensate. In stud- 
ies of the condensate it is important to know the low tem- 
perature elastic scattering cross sections, which control 
the thermalisation of the atoms. The spin-change cross 
sections are also important. Atoms interacting through 
the triplet state of NaRb remain in a magnetic trap but 
those in the singlet state leave. Thus spin-change is an 
important trap loss mechanism. We present calculations 
of the scattering lengths and volumes and the effective 
ranges that suffice to quantify low energy elastic scatter- 
ing, and of the spin-change cross section, which quantifies 
trap loss. We note that in a full treatment of stability 
of the condensate the hyperfme structure must be ac- 
counted for 0, 0] but in such a treatment the scattering 
data that we present are needed. 

The two main processes that dictate the properties of 
a trapped cold atomic gas in thermal equilibrium, are 



'Electronic address: mjj@dcs.gla.ac.uk 



the elastic collision and the spin change collision. The 
clastic collision process is essentially of s-wave nature if 
the temperature is sufficiently low. The spin change col- 
lision is an inelastic process during which the spins up 
and down (s and s') of two scattering atoms, A\ and A 2 , 
are exchanged: 



A 1 (s) + A 2 (s')~ *A 1 (s')+A 2 (s) 

In addition to the stability criterion that has to be sat- 
isfied to obtain a BEC, it is also important for practical 
purposes to ensure that the system allows efficient evap- 
orative cooling. This is possible if the ultracold s-wave 
clastic collisions dominate the inelastic spin change col- 
lisions, i.e. if the related cross sections satisfy a c \ 3> cr sc . 
For a given pair of atoms of different species in the sin- 
glet or triplet state (s/t), the low energy elastic scattering 
cross section is proportional to the square of the s-wave 
scattering length, and the spin change cross sec- 

tion, to the square of the difference of the two scattering 
lengths, oi^ — as . It is of interest to estimate an up- 
per bound to the temperature range in which only s-wave 
scattering is significant and to investigate the dependence 
of the spin-change cross section on the atomic masses and 
the temperature. 

Previously we demonstrated the utility of the vari- 
able phase method to circumvent numerical technical dif- 
ficulties and to accelerate convergence of the accumu- 
lated scattering parameters [1, Q. However the relia- 
bility of computed data also depends on the accuracy 
of the molecular potentials; it is possible that, within its 
error bounds, a molecular potential supports a near zero- 
energy bound state which makes the computed scattering 
length extremely sensitive to any change in the potential. 
Cold collisions of Na and Rb atoms are interesting in pho- 
toassociation and trap loss experiments and we evaluated 
scattering lengths and volumes, and numbers of bound 
states computed using various potentials @. We inves- 
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tigated the mass dependence of our results, which were 
compared to the data published by Weiss et al. It 
has been suggested by Cote et al. [ll[ that our calcu- 
lations were tentative because the potentials need to be 
specified more precisely. The new potentials for NaRb of 
Pashov et al. [121 ] represent an improvement over those 
that were available at the time of our calculations but, 
contrary to their suggestion 11211. the confusion of units 
in the work of Weiss et al. [H| docs not compromise 
the reliability of our earlier calculations since we did not 
use but merely quoted their results, which were corrected 
later 

For the present work, we constructed molecular poten- 
tials using the improved short range data of Pashov et al. 
fl2| . The exchange potential is given by the expression 
of Smirnov and Chibisov [lij in the formulation used by 
Weiss et al. ■ Smirnov and Chibisov [3] showed that 
the exchange potential between two heteronuclear atoms 
in their ground ns 2 S states can be approximated by the 
series (which having a single term is exact for homonu- 
clear atoms) 



van der Waals coefficients in proportion by amounts of 
order of the inverse of the reduced mass, but in these cal- 
culations these corrections would be smaller than 0.01% 
and therefore we do not include them. 

In Section 2, we give an overview of the theory and dis- 
cuss the semiclassical approximation to the p-wave scat- 
tering volume. We show how we produce estimates of 
the temperature range in which only s-wave scattering is 
significant from a simplified analysis of the s-wave and 
p-wave cross sections; and we derive an expression for 
the low energy spin-change cross section. In Section 3, 
we present and discuss numerical data. 



II. THEORY 

A. Low energy scattering 

The accumulated s-wave phase, 5k{R), for wavenumbcr 
k of relative motion, satisfies the differential equation in 
the nuclear separation R of the Riccati type [23|] : 



N 



V E (R) = -J2 



J n R! 



i? 7 exp(-/3i?), (1) 



where Ve(R) hartrec is the exchange potential at R 
bohr, cn 2 /2 hartree is the ionisation energy of an atom, 
(3 = a + a' and 7 = 2/ a + 2/ a' + l/(a + a') - 1; J„ is 
calculated from the asymptotic behaviours of the wave- 
functions of the valence electrons of the atoms and it 
also depends on quadratures over functions with branch 
points. We evaluated the values of J„, shown in Table 
HI for alkali metal atom pairs using the calculations of 
Marinescu and Dalgarno [l6[ for the asymptotic wave- 
functions of the valence electrons. The definition of J n 
is half that used by Smirnov and Chibisov [H| ■ We used 
experimental ionisation energies [13] ■ The values of Jn , 
Ji and J2 for NaRb agree with those of Weiss et al. [Hi • 
We note that Zemke and Stwalley [HI proposed a rep- 
resentation of the exchange interaction equivalent to a 
single term of equation (1), —CR a exp(— (3R) with a and 
(3 to be determined. For consistency in comparison with 
the calculation of Weiss et al. we used their formulation 
i.e. equation (1). 

Long range molecular potentials of alkali-metal atoms 
are usually well characterised by three van der Waals 
dispersion coefficients Cg, Cs and Cio fl8lh we took the 
values computed by Derevianko et al. [l9l. l20l | and we 
added higher order inverse powers where we used the co- 
efficients beyond Cio computed by Mitroy and Bromley 
plj . The short range potentials of Pashov et al. fl2j . 
being experimentally based, account for adiabatic, dia- 
batic and relativistic effects. The long range dispersion 
potentials do not account for such effects but we can 
consider adiabatic corrections. Dalgarno and McCarroll 
p2| showed that allowance for the adiabatic coupling be- 
tween the nuclear and the electronic motion changes the 



dS k (R) 
dR 



= -fc _1 y(i?) sin 2 [kR + S k {R)} 



(2) 



where V(R) is 2^,/h 2 times the interaction potential and 
\i is the reduced mass of the dimer. The phase shift is 
the limit of the accumulated phase at infinite separation 
and we write it as 8k = 8^(00); it suffers no mod[7r] am- 
biguity. The number of s-wave bound states supported 
by the potential is given by Levinson's theorem pH, [24| 
as 7riVb = lim fc ^o 4- 

In effective range theory (2f| HH, [27| the low energy 
expansion of the tangent of the s-wave scattering phase 
shift yields the equations satisfied by the accumulated 
scattering length, a s (R), and a volume, b s (R), that is the 
coefficient that appears to third order in the expansion 
of tani5fc(i?) for small values of k: 



das(R) = [R-a s (R)] 2 V(R) 



dR 
db s (R) 
dR 



(3) 



[R-a s (R)Y 



-2[R-a s (R)}b s (R)\V(R!}i) 



The accumulated effective range, R c ff(R), is related to 
the accumulated scattering length and the volume by: 



Rctf(R) = 



2a s (R) 2b s (R) 



[a s (R)Y 



(5) 



Eq. (HJ), satisfied by b s (R), is coupled to Eq. and has a 
closed form solution [23, [28| but the solution of Eq. (J4| is 
difficult to compute because of the divergence problems 
that reflect the presence of poles in the accumulated scat- 
tering length a s (R). We have discussed this in detail in 
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Ref. [28|, where an analytical first order long range cor- 
rection to the effective range is given. In Refs. [a Q we 
showed how upper and lower bounds to the scattering 
length and scattering volume can be obtained and gave 
expressions of improved approximations to the scattering 
parameters as linear combinations of those bounds. 

In the following we mean by the radial wave func- 
tion, R times the actual radial wave function and, by 
the Schrodingcr equation, the equation that is satisfied 
by this radial wave function (that contains no first deriva- 
tive). In our previous analysis [§] we related ^-wave 
scattering to s-wave scattering using a relation given by 
Calogero [23|; here we present a simplified version for p- 
wave scattering. Making the substitution Z = i? 3 /3 and 
writing the zero energy radial wave function as R~ 1 (j){Z), 
we obtain, from the p-wave Schrodingcr equation, an s- 
wave Schrodingcr equation that is satisfied by 4>{Z) but 
with the potential divided by (dZ/dR) 2 . The asymptotic 
solution, suitably normalized, is Z— a where a is a volume 
analogous to the scattering length. Application of vari- 
able phase theory with 4>{Z) = Z — a(Z), where a(Z) is 
analogous to the accumulated scattering length, yields an 
equation like Eq. ©. However the function R^ 1 (Z — a) 
is the same as the function R 2 /3 — a/R, which is the suit- 
ably normalised asymptotic p-wave function (29[ . Hence 
the s-wave scattering quantity a for the potential divided 
by (dZ/dR) 2 is the p-wave scattering volume. On writ- 
ing the equation in Z that is equivalent to equation ([3]) 
we find the following equation that is satisfied by the 
accumulated p-wave scattering volume a p (R): 



da p (R) 
dR 



fR 3 



- ap(R) 



R~ 2 V(R). 



(6) 



as R~ n . Below we compare our quantal results for the 
scattering volume with those predicted by Dickinson's 
formula. 



B. Temperature range 

Knowledge, even approximate, of an upper bound for 
the temperature range, T, in which s-wave scattering 
dominates the elastic collisions in the trapped gas of ul- 
tracold atoms, is useful. Here, we show how estimates 
can be found from a simple analysis of the s-wave and 
p-wave scattering cross sections. For small values of the 
wavenumber k the s-wave scattering contributes 
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, 2 a 2 



k R cS a s 



(8) 



to the elastic cross section. If the gas is sufficiently cold 
that fc|a s | <§; 1 and fc 2 i? c ff|a s | <C 1, where k now denotes 
the thermally averaged wavenumber, then the cross sec- 
tion is not influenced by the effective range and is also 
well approximated by <7 e i,s ~ Aira 2 . If the second inequal- 
ity, k 2 R c ff \a s \ <C 1, is satisfied then so is the first (unless 
the effective range is pathologically large). In this case, 
the temperature T satisfies T <C h 2 /2/ikBR c s \a s \- In esti- 
mating an upper bound to the temperature, let us require 
that the value of k be smaller than 1% of 1/ -J R c g \a s \. 
The contribution of p-wave scattering to the elastic cross 
section is given by a e \ tP = 

ten as c e i )P = Ana 2 x 3fc 



127rfc 4 a 2 , which can be rewrit- 

- - 2 

I . If k satisfies the 



(a p /< 



above criteria, which it does if T -C h 2 /2/j,kBR c s |a s | , then 
<T e i iP /er c i,s <C 1 and p-wave scattering can be neglected. 



Recently Dickinson [36( derived a semiclassical formula 
for the scattering volume. We can derive his formula 
directly from the relation 



poo f oo i r> 

/ ^V{R) dR = / ^V{Z) — dZ (7) 
Jro Jz dZ 

where V(Rq) = and Z a = i?o 3 /3, and the semiclas- 
sical s-wave formula of Gribakin and Flambaum [37| . 
With Z = R 3 /3, the leading term of the long range po- 
tential is -(C 6 /27s/3)/Z 10 / 3 where C 6 is the van der 
Waals dispersion coefficient (in terms of R the leading 
long range term is —Ce/R 6 ). Substituting a = Cq/27-\/3 
and n = 10/3 into the formula of Gribakin and Flam- 
baum we obtain the semiclassical formula of Dickinson. 
The Schrodingcr equation for ^-wave scattering can be 
reduced to an s-wave scattering equation by the transfor- 
mation Z = R 2l+1 /(2l + 1) [alia]. This transformation 
and the formula of Gribakin and Flambaum can be used 
in the manner described above to find a generalisation of 
the semiclassical formula for the quantity analogous to 
the scattering length, provided it exists, for higher angu- 
lar momenta and potentials that behave asymptotically 



C. Spin change cross section 

The spin change cross section is [30} : 



l)sin 2 



°k,l 



k.l 



(9) 



where the 6^] and denote the phase shifts of the 
partial wave with angular momentum I, in the singlet 
and triplet states respectively (note that for the s-wave 
case, in Eq. @, neither the ^-dependence of the phase 
shift nor the spin configuration, were made explicit). 

In the ultra- low energy limit {k ~ 0), only the partial 
waves with angular momenta I = and I = 1 contribute 
significantly to the scattered wave function, the contri- 
butions from higher angular momenta being negligible at 
very low temperature [27l . l3ll . l32l . l33l |34| . The low energy 
expansion of the s-wave phase shifts, in the singlet/triplet 
states, is: 
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(10) 
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and 



b,l 



(11) 



for the p-wave phase shifts. and are the 

numbers of s-wave and p-wave bound states. The volume 
parameter bi s ^^ is related to the s-wave scattering length 
and effective range through equation ([5]). Note that in 
the limit k — > 0, Levinson's theorem is recovered from 
Eqs. (dSJ) and {TTJ). 

Now, inserting Eqs. (fTO)) and (fTTj) into Eq. ([9]), we ob- 
tain the following expression for the spin-change cross 
section, in the low energy limit, up to second order in 
relative motion energy (i.e. up to fc 4 ): 



a sc , = 47r 



bP - b^ 



4 8) - 4 ] 



bP - b^ 



k 2 

(12) 



The first term of the right hand side of equation (fl~2| , 

2 



4tt 



(t) 



, represents the zeroth order ap- 
proximation of the spin change cross section in the low 
energy limit, and it is temperature independent. The 
spin-change cross section, a sc , is studied as a function of 
isotope mass and cold atom system temperature, in the 
next section. 



TABLE I: a and J n (atomic units). 



10 4 x J n 

n Li Na K Rb Cs 



Li 



0.62951 
1 

2 
3 
4 



172.85 



K 0.56483 
1 
2 



Rb 0.55409 
1 

2 



152.22 
1.4498 
21.500 



89.522 
3.9634 
12.475 
1.3507 
4.0052 



80.350 63.424 

4.2131 4.2858 

11.184 8.8276 

1.4304 1.4460 

3.5883 2.8334 



Na 0.61459 134.13 79.050 70.986 56.082 

1 2.7489 3.0487 3.2593 

2 10.900 9.7686 7.7056 

3 0.93017 1.0275 1.0911 

4 3.4689 3.1046 2.4470 



46.942 42.226 33.466 
0.34607 0.78409 
5.6200 4.4215 



37.999 30.138 
0.45928 
3.9529 



Cs 0.53497 



23.936 



III. RESULTS AND DISCUSSION 

First, in Table HI we give the coeffients a and J n , 
appearing in Eq. (JTJ) , computed for alkali dimers. The 
masses that we used in our calculations of the scattering 
data are 22.98976967 for 23 Na, 84.9117893 for 85 Rb and 
86.9091835 for 87 Rb (in atomic mass units) [35[. 

The scattering data given in Table [IT] arc calculated 
from the potentials in Ref. (T^]. A useful check on the 
reliability of our numerical procedures is the agreement 
of our s-wave scattering lengths with those of Pashov et 
al. [HJ for both isotopomers in the singlet and triplet 
states; new data that we provide are the p-wave scat- 
tering volumes, s-wave effective ranges and numbers of 
bound states. Table ILT1 shows that the values of the scat- 
tering volumes are close to those predicted by Dickinson's 
relation between the scattering length and the scattering 
volume (Eq. (21) of Ref. (36|) except for the 87 Rb isotope 
where the quantal and semiclassical values, a p and ap C , 
are not so close themselves. Note that both values remain 
large in comparison to the other volumes, which, together 
with a change of sign, possibly indicates the presence of 
a p-wave resonance for the 23 Na - 87 Rb cold collision in 
the singlet state. The presence of a resonance, which is 
further investigated below, explains the sensitivity of the 
scattering data to the semiclassical approximation. 



Despite the sensitivity of the computed scattering data 
to minor changes in the molecular potentials, we find that 
influence of the van der Waals terms with coefficients Cn 
to Ci6 is negligible: they contribute less than 1% to the 
computed scattering parameters. This is consistent with 
the conclusion of Mitroy and Bromley [2l[ that, because 
of cancellation, the contribution of the dispersion terms 
—C n /R n might be neglected when n > 10 if the accuracy 
of the van der Waals coefficients is of the order of 1%. 
The temperatures below which only s-wave scattering is 
significant are all smaller than 100 nK; thus p-wave scat- 
tering can be neglected only in extremely cold Na and Rb 
gases, particularly for gases containing the 85 Rb isotope. 

We turn now to the study of the spin change cross sec- 
tion, whose main contribution comes from cr® c . We ob- 
tain: ct° c = 1.25 x 10 6 bohr 2 and cr s ° c = 1.91 x 10 4 bohr 2 
for the 23 Na 85 Rb and 23 Na 87 Rb isotopes respectively 
Considering the values of the elastic cross sections given 
in Table [TH the first important observation is that the 
criterion a c \ 3> er sc is not fully satisfied for both iso- 
topomers: stability of the dual species 23 Na 85 Rb and 
23 Na 87 Rb condensate against the inelastic spin change 
process can hardly be achieved. It is also interesting to 
note that the spin change cross section is also very sen- 
sitive to a variation of the atomic mass since it is two 



TABLE II: S-wave scattering lengths a a , elastic scattering 
cross sections <t c i,s, p-wave volumes a p (quantal) and a^ 
(semiclassical), effective ranges R c s, s-wave and p-wave num- 
bers of bound states Nb,o and Nh,i, and temperatures T. 



A 1 ^ 



3v+ 



a 3 E 





23 Na 85 Rb 


23 Na 87 Rb 


23 Na 85 Rb 


23 Na 87 Rb 


a s (bohr) 


396 


109 


81 


70 


rr ■ (hnhr 2 ) 


1 Q7x10 6 


1.49x 10 s 


a o % v l n 4 


is i fix i n 4 

U.IDA ±\J 


a p (bohr 3 ) 


-4.30 x 10 5 


7.0 x 10 6 


3.10 x 10 5 


1.30 x 10 5 


al c (bohr 3 ) 


-4.25 x 10 5 


1.38 x 10 7 


3.15 x 10 5 


1.31 x 10 5 


ReS (bohr) 


31 


58 


84 


113 


A b , 


83 


83 


23 


23 


N hA 


82 


83 


23 


23 


T(nK) 


53 


75 


70 


60 



orders of magnitude larger for the lighter system than 
it is for the heavier one. This is coherent with the fact 
that because of the lighter mass of the atom 85 Rb, the 
larger spatial extension of its wavefunction enhances the 
electron spin exchange during the scattering act. 

On the left panel of Fig. [TJ the variation of a sc (T) — <7g C 
is shown for temperatures up to the /iK domain. It is ob- 
vious that in this temperature range the first and second 
order corrections contribute only a little to the total spin 
change cross section. These corrections remain larger in 
the case of the lighter isotopomer for all temperatures in 
the range studied. However, if we define a standardized 
spin change cross section as ct sc (T) = [a sc (T) — <j° c ] /<7° c , 
shown on the right panel of Fig. [1] we find that for tem- 
peratures close to 60 K to 70 K, the behavior of the cor- 
rections changes: the amplitude of a sc (T) for 23 Na 87 Rb 
becomes larger than it is for 23 Na 85 Rb. This is due to the 
fact that the last term of the right hand side of Eq. (fl2)) . 
proportional to the square of the difference of the p- 
wave volumes, becomes relatively important (compared 
to the other first- and second-order correction terms) as 
the temperature increases; specially considering that the 
magnitude of the p-wave volume of 23 Na 87 Rb in the sin- 
glet state, is particularly large compared to the other 
volumes. Note that this finding is coherent with the 
simple analysis of the s-wave and p-wave elastic cross 
sections we presented above: p-wave scattering becomes 
non-negligible in this temperature range. 

Furthermore, in the case of the singlet state charac- 
terised by the Af 1 S + molecular potential, the signs of 
the scattering volumes for the two isotopomers differ. To 
interpret this result wc first compare the numbers of p- 
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FIG. 1: Variations of cr BC (T) - cr s c (bohr 2 ) and a sc (T) with 
temperature T (K). See text. 




FIG. 2: Variation of p-wave volume a p (10 8 bohr 3 ) with re- 
duced mass. See text. 



wave bound states of the two isotopomers. These can 
be obtained by applying Levinson's theorem in the same 
fashion as for s- waves. The accumulated p-wave phase 
shift, <5fc,i(i?), satisfies the equation [23[: 



dS kA (R) = V(R) 
dR k 



cos 5k,i{R)j(kR) 



wa.5 k>1 (R)h(kR)Y 
(13) 

where j and n denote the Riccati-Bessel functions of or- 
der 1: j(z) + ih{z) = (iz~ l — l)cxp(— iz). The number 
of p-wave bound states supported by the A" 1 E + poten- 
tial is changed from 82 for 23 Na 85 Rb to 83 for 23 Na 87 Rb, 
demonstrating that with the increased reduced mass of 
the heavier isotope the interaction potential supports an 
extra bound level. In the case of the triplet state char- 
acterised by the a 3 E + molecular potential, the numbers 
of p-wave bound states are the same, 23, for both iso- 
topomers. The much larger value for 23 Na 87 Rb suggests 
the possibility of a near zero energy bound state. 

To find evidence for this we studied the dependence 
of the p-wave volume on the reduced mass of the NaRb 
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molecule. We have defined a reduced mass, fjf, which 
varies between that of the lighter isotopomer, fxi, and 
that of the heavier one, [i^: 

n' = m + \(fi h - m), (14) 

where A is in the range [0,1]. The resonance in the varia- 
tion of the scattering volume with reduced mass, shown 
in Fig. [21 suggests that the potential just supports a p- 
wave bound state for the heavier isotope. At resonance 
the parameter A ~ 0.9625. The s-wave scattering length 
at the corresponding mass is 116 bohr and, as predicted 
by Dickinson [3(| , is approximately twice the mean scat- 
tering length defined by Gribakin and Flambaum [37l |. 
which is 55 bohr. Table lU shows that the values of 
the scattering volumes are close to those predicted by 
Dickinson's relation between the scattering length and 
the scattering volume [36[ except near the resonance. In 
Ref. [36[ , the quantal and semiclassical values are in good 
agreement because the molecular potentials used for the 
calculations, do not support an additional p-wave bound 
state as does the improved potential of Ref. [l2j ■ 

IV. CONCLUSION 

We calculated parameters that describe low energy col- 
lisions of Na and Rb atoms. We confirmed that the long 



range part of molecular potentials is sufficiently well char- 
acterised by the set of three van der Waals dispersion 
coefficients Cg, Cs and Cio to calculate the scattering 
data. We estimated temperatures below which only pure 
s-wave contributes to the scattering cross sections and 
showed that the inelastic spin change collisions dominate 
the s-wave elastic scattering in these systems. We found 
evidence for the existence of a near zero energy p-wave 
bound state. This is in accord with the semiclassical 
study made by Dickinson [36| . 
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